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Abstract 

A unit volume zero intelligence (ZI) model is defined and the distribution of its LI 
process is recursively described. Further, a generalized ZI (GZI) model allowing non-unit 
market orders, shifts of quotes and general in-spread events is proposed and a formula 
for the conditional distribution of its quotes is given, together with a formula for price 
impact. For both the models, MLE estimators are formulated and shown to be consistent 
and asymptotically normal. Consequently, the estimators are applied to data of six US 
stocks from nine electronic markets. It is found that more complex variants of the models, 
despite being significant, do not give considerably better predictions than their simple 
versions with constant intensities. 
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1 Introduction 


With the recent wide expansion of trading according to continuous double auction (CDA), the 
importance of mathematical modelling of this trading mechanism grows. 

A number of models y>f CDA exist which assume a rational behaviour of agents involved 
Parlour and Seppi f 2008f) and the references therein); these models are, however, 


fsee, e.g., 

dependent on many arbitrary assumptions and do not give much better e mpirical r esu lts than 


see 


Gode and Sundeij ( 1993 1 for 


models assuming a purely random behaviour of agents 
discussion). Thus, as an alternative to the “rational” approach, zero intelligence (ZI) models 
of the C DA s t arted to b e stu di ed: out of a number of similar m o dels of that kin d, let us name 
Stieler ( 1964h : Maslov f 200(1 ; Challet. and Stinchcombel ( 2001 1; Luckockl ( 20031) ; Smith et al. 


( 2003 1: Mike and Farmer ( 2008 1 or Cont et al. (j2010h . All these models assume unit order sizes, 
Poisson arrivals of market and limit orders and locally constant cancellation rates depending 
on the distance to the qu otes !)] For a survey of the ZI m odels and their characteristics, see 
Chakraborti et al. ( 201 lah and Chakraborti et al. ( 2011b), e specially Section If of the latter 
paper. For the recent developments, see lAbergel et al. (201 ill. 

Although, according to their advocates, the ZI models are able to mimic many stylised 
pirical facts such as fat tails or non-Gaussianity ( Slaninal ( 200ll 2008 1: Snrid ( 2008 1: Cont 


em 


(2011)), no rigorous statistical evidence in that respect has been presented yet due to intractabil¬ 
ity of these models. In ISmfdl (2012), a conditional distribution given LI data is described; 
however, the model considered in this paper is too general to ensure consistency of statistical 
estimators. 

The present paper simplifies the approach of Smfdl ( 2012tl by assuming only a finite number 
of prices so that it is possible to construct consistent asymptotically normal estimators of the 
models’ parameters. 

The exposition is started by introducing a sufficiently general set ting cove ring a wide range 
of existing zero in telligence (ZI) models, e .g., an early wo rk of IStiglerl ( 196411 . the model by 
Cont et al. ( 2010h . a discretised version of Luckock ( 2003h and a slightly modified version of 
Sm ith et al. (12003 b After demonstrating ergodicity of the covered models, a conditional density 
of jumps of the LI process (i.e., the bid, the ask, and the corresponding offered volumes) given 
its history is formulated J§ 

Further, in order to treat the most obvious discrepancies between the ZI models and reality, 
a generalised (GZI) setting is proposed, which allows shifts of the quotes, non-unit market orders 
and general distributions of inside-the-spread events. A formula for a conditional density of the 
out-of-spread jumps of the quotes, later used in estimating the in-book parameters, is given. 

Subsequently, Maximum Likelihood estimators are formulated and several variants of both 
the ZI and GZI models are estimated by means of LI data for six stocks on nine US electronic 
markets. Further, for each stock-market pair, the variants of the models are compared by their 


1 The only exceptions are Maslov ( 2000l) and IChallet and Stinchcombe (| 2001 ) with a discrete time (which 
could, however, be regarded as Poisson events’ time) and Mike and Farmer! (|2008h . in which the cancellation 
rate also depends on the order books’ size and imbalance 

2 Although order book (L2) data might seem more natural to be used for the estimation, the better availability 
and lower price of LI data speak for using them. This approach, in addition, allows us to avoid a problem of 
hidden limit orders, which are invisible in the books but affect sizes of quotes’ jumps. 
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ability to forecast magnitudes of the quotes’ jumps. It is found that more complicated variants 
of the models do not bring significant improvements in comparison with the simple variant with 
constant intensities, which itself, however, nearly always performs significantly better than a 
naive prediction of the jumps. It is also found that GZI variants of the models do not perform 
significantly better than their ZI counterparts. 

This paper is organised as follows: First, the ZI setting is defined (Section [2]) and the 
distribution of the LI process given its history is derived. Consequently, the GZI setting is 
introduced and a formula for the density of the out-of-spread jumps and volume changes is 
given (Section [3]). Finally, the empirical evidence is discussed (Section Hj) and the paper is 
concluded (Section [5]). The Appendix includes a proof of Theorem 12.21 (Appendix [All . a proof 
of asymptotic properties of the MLE estimator (Appendix [B]), and detailed results on the 
estimation (Appendix £Tj) . 


2 Zero Intelligence Model 

2.1 Definition 

Consider a general discrete-price continuous-time zero intelligence model with unit order sizes 
described by a pure jump type process 


— (^4* 5 B t )i t > 0, 

where 

A t = (A],Al...,A?)eN n 0 , 

and 

B, =(Bi,Bl...,B?)eNS, 

are the sell limit order book and buy limit order book, respectively; here, n is a number of 
possible prices (ticks) and, for any p, Af (T>f) stands for the number of the sell (buy) limit 
orders with price p waiting at time t. Further, denote by 

a t := inf {p : A t {p) > 0} V 0, b t : = sup {p : B t (p ) > 0} A n + 1, 

the (best) ask and bid respectively. The list of possible events causing jumps of S together 
with the notation for their intensities is given in the following table: 
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code - e 

intensity - I t (e ) 

description 

BMO 

9 t = 9(at, bt) 

An arrival of a buy market order, causing 
A at to decrease by one (if the sell limit order 
book is empty then the arrival of the market 
order has no effect). 

SLO(p) 

Kt, P = K(at,b t ,p) 

An arrival of a sell limit order with limit price 

p > b t causing an increase of A p by one. 

CA(p) 

Pt, v = A p t p(a t ,b t ,p ) 

A cancellation of a pending sell limit order 
with a limit price p causing a decrease of A p 
by one. 

SMO 

dt = d(a t , b t ) 

An arrival of a sell market order, causing B bt 
to decrease by one (if the buy limit order 
book is empty then the arrival of a market 
order has no effect). 

BLO(p) 

A t ,p = A (a t ,bt,p) 

An arrival of a buy limit order with limit 
price p < a t causing an increase of B p by 

one. 

CB(p) 

cr t ,p = Bf<j(a t ,bt,p) 

A cancellation of a pending buy limit order 
with a limit price p causing a decrease of B p 
by one. 


Here, all 9, k, p, d, A, a are measurable functions. 

It is assumed that all the flows of the market orders, and the flows of limit orders as well as 
their cancellations are mutually independent in the sense that the conditional distribution of 
the relative time of the first event following t G given the history of S up to t is exponential 
with parameter 


n at — 1 n bt 

At — 9 t + 'd t + K Pt t + \>.t + Bfcr P! t 

p=bt +1 p =0 p=a t p =0 


( 1 ) 


and that 


the probability that the next event will be of type e equals to I t (e)/A t (2) 

where It(e) is the intensity of event e at t. It is obvious that S is then a Markov chain in a 
countable state space. 

Finally, it is assumed that ao,&o are deterministic and 

Hq 0+1 , Aq 0+2 , ... are Poisson with parameters i ao +i , t ao + 2 , • • • , 

mutually independent and independent of B 0 (3) 

where l, > 0 are constants, and that a symmetric assumption holds for B. 


2.2 Relation to Existing Models 


The following table shows how some of the m odels mentioned in the Introduction comply with 
our setting. When speaking about Luckock (2003]), we have its discretised version (see Smfd 
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(2012), Sec 3.3) in mind. When speaking about Smith et al. (120031 ). we are considering its 
bounded version (i.e., contrary to Smith et al (120031 ) we assume zero arrival intensities for 
prices smaller than one and greater than n ). 


model 

e 

■d 

K 

P 

A 

a 

Luckock (2003) 

K(b t ) 

1 - L(a t - 1) 

K L p 

0 

A^ 

0 

Smith et al. (2003) 

6 s 

6 s 

K S 

p s 

K S 

P s 

Cont et al. (2010) 

9 C 

9 C 

K c (p - b t ) 

•+J 

-o 

1 

3 

O 

n c (a t - p ) 

P c (a t — p) 


Here, n l p = K(p) — K(p — 1), \ l p = L(p) — L(p — 1) where K , L are (continuous) cumulative 
distribution functions, n c , p c are measurable functions and the rest of the symbols are constants. 

Some of t he mode l s from th e Introduction were not mention ed in the table: We did not 
include either iMaslovi (1200011 or IChallet and Stinchcombel (1200111 because they both conside r 
discrete time and are very s imilar to ISmith et al. (120031) with p = a = 0, ICont et al. (120101) . 
respectively. The model by Mike and Farmerl (1200811 was not included because of its compli- 
cate d cancellation sub-model and because, apart from the cancellations, it is similar to that of 


Cont et al, (2010). Finally, we did not include ISti glerl (19641 1 because it is a special version of 
Luckock ( 2003h (with Ii(x ) = L(x) = x). 


2.3 Distribution 

Let us start with a result which will guarantee that sampling from the model will give enough 
information to a statistician. 


Proposition 2.1. If />(•) > 0, cr(«) > 0, $(•) > 0 and $(•) > 0. then S t is ergodic. 


Proof. The Proposition may be proved analogously to Cont et al. ( 201 0). Proposition 2, where 
the ergodicity is verified by finding a Markov chain in No dominating the total number of orders 
with ajrecurrent zero state. In particular, it suffices to replace the definition of A and /i, from 


Co nt et al (2 0101 ) by 


A=p [max n(a, b,p) + max A (a, b,p)\, 

* ^ ci,b,p a,b,p 

P= 1 


and 


Pi = min 9(a,b) + mind(a,6) + i min[min(p(a, b,p), a(a, b,p)\. 

a,b a,b a,b,p 


□ 


Our next goal is to derive a recursive analytic expression of the distribution of 

x T = (a T , h T , q T , ry) T > 0 , q T = A a T T , ry = B b r T , r > 0 

(the LI process). To do so, let us denote by • • • the jumps of x and, in order to avoid 
frequent double indexing, let us write instead of y ti and instead of y ti - for any process 
y. For each i, we want to evaluate 


P[(Atj,X[i]) G •|f[ i -i]], £[k\ = (x[k], Atk,x[ k -i],...,x 0 ), k> 0, 
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starting with 


P[(Ati,a w , ?[,-]) G •!%-!]. • 


To this end, note that 

• a jumps down if and only if a limit order arrives into the spread, in which case a limit 
price of the new order becomes a new value of a. 

• a jumps up if and only if the offered volume of the ask decreases to zero clue to either a 
market order arrival or a cancellation, in which case the closest occupied tick becomes a 
new value of a. 


Formally, 

(^[i-l] + I; 1) e i — SLO(f>[j_!] + 1) 

(^[i-i] + 2,1) ei = SLO(6[i_i] + 2) 

(®[i—i] 1) 1) c* SLO(a[j_i] 1) 

( a [i-lp Q[i-1] + 1) e i = SLO(a[j_i]) 

(«[*—i], Q[i-i] — 1) e i £ (BMO,CA(a[j_i])},g[j_i] > 1 
G (BMO,CA(a[j_i])},g[j_i] = 1 
(a.p_i], gp_i]) otherwise 

di = inf{7r : 7r > ap_i], Ay > 0} A n + 1 

where e, is the type of an event happening at A; and where A n+1 = 0 by definition. Thus, to 
determine the conditional distribution of (Atj, ayq, q[{\), it suffices to know a joint distribution 
of (Atj, e*, A[j_]) which is described by the following Theorem in three steps: 

Theorem 2.2. (%) For any r and e, 



where 


and, for any k, 


P[Atj > T,e t = e|%_ij] = exp{-7j_ir}7rj(e), 

^p\i— 1]> ^ SLO (jp) , byi— j] p A 0[j_i], 

%- 1 ], e = BMO, 

w,(e) = -y~\ x h"- 1 '' e = CA (“l‘->l). 

A p> [i_i], e = BLO(p), 6[i_ij < p < ap—ij, 

%_!], e = SMO, 

k o-[*-i] e = CB(6[j_ij), 

a [fc] a [fc] — i 

Tfc ^[fc] T T P[fc] T T ^ ^ A ^ ^ Ap,[fe]- 

p=6[fc]+ 1 P=b[k] 
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(it) Denote by S c the Dirac measure concentrated in c (i.e., the distribution of constant c) and 
write o for convolution (i.e., summation of two independent variables). For any i and p, 


!-]!£[*-] 

where 
f[H = (e i; 


0 on [1 < p < a[i- 1 ]] 

on [P = a [i- 1 ]] 

Bi (v P!i , zu P}i ) o Po (e P)i + t P)i ) on [a[j_i] < p < n] 


^p,i 


= f A [ki(p)-] ^) >0 > 

0 ki(p) = 0, 


= j n;=U) fa h{p)<i 

1 ki(p) = i 




Cp,i ~ 


i— 1 

'y ] 0p,j(i _ $ 

j=h(p) 


i— 1 

VO) dp,m, 

m=j-\-1 


Lp,i 


k t {p) > 0, 

h> IL=o dp,j hip) = 0, 


hip) = minjfc > 0, a^j < p] V i 


and, for any k, 


Pp,k — 


^p,[k\ 

Pp,[k\ 


d P ,k = exp{-p Pt[k] At k+1 } 


(m) _j are mutually conditionally independent given £[,_] 

Proof. See Appendix 1X1 


□ 


Before proceeding further, let us illustrate point (ii) of the Theorem, which is somewhat opaque, 
by an example. 

Exapmle 2.3. Consider model by Smith et al. (20Oil) with i. = 0 (implying t #3 = Oj. Let 

1 = 3, a 0 = 2, q 0 = 1, ap] = 3, qp = 2, ap] = 1, q[ 2 ] = 3 (see FigureUlfor an illustration). 

Let us start with p = 2: as no orders could be present in tick 2 before t 2 (the ask was above 

2 that time), all the in-book orders present in the tick f 3 are those having arrived between t 2 
and t 3 and not being cancelled. According to Appendix\A J number of those orders is Po(e; ;3 ) (to 
check it, note that h(3) = 2 so v l2 = A^ 2 _j = 0 and e,; )3 = ^[1 — exp{—pAt 3 }] which, according 
to Appendix\A 1 is actually the distribution of the immigration and death process of length A t 3 ). 

For p = 3, the situation is the same with the difference that number of initial orders having 

survived until t 3 , which is distributed according to Bi(2, exp{—pAt 3 }) = Bi(v4^_j, <5*, 3 ) ? has to 

be added to Af ,. 

[H 

Further, as the orders present in tick f at t 3 are exactly those having arrived from the start 
and not being cancelled, A|_j|£[j_] ~ Po ^[1 — exp{—pf 3 }j - it is easy to check that (ii) gives 
the same result. 

Finally, because fcj(l) = 3, we have = 0, = 1 and u\ = 3 so, according to (ii), 

A\_ ,~Bi(3,l) = tf 3 . 


From the Theorem and (HI) , we get that 
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Figure 1: Illustration of Example 12.31 


Corollary 2.4. (i) Conditional density of (e*, Atj|£[j_!]) is given by 

fi{ T i e |£[i-i]) = lff\ exp{—7i_ir}7Tj(e), 
(ii) conditional density of (ap], Q'p]) ] is given by 


gi(a,q\£[i-]) 

1 [q = + 1] on [6[j_i] < a < a [*_!], e* = SLO(a)] 

1 [q = q\i- 1 ] - 1] on [a = a[i_ 1 ], Cj e {CA(o[i_i]), BMO}, > 1] 

_ a— 1 

w a ,,,i El (! " on [a.p_i] < a, e* G {CA(a [t _i ] ),BMO},g [i _ 1] = 1] 

lr = a [i-l]+ 1 

l[g = a = otherwise , 


where 1 denotes indicator function and 


^p,q,i 


P[Bin7 Pli ) o Po(e Pii + t Pii ) = q\ < p < n, 

l[q = 0] a = n + 1. 


(Hi) conditional density ofia^q^e^/Sif) given £[i-i] is given by 

hi(a , q, e, = gi(a, q\e, t, £[i_i])/;(e, [*_!]). 

Remark 1. As the definitions of B and A are symmetric, formulas for the (conditional) distri¬ 
bution of (b, r) are symmetric to those for (a,q). Moreover, as changes of (a,q) and (b,r) are 
mutually exclusive almost sure, i.e., 

[(a, q) changes at t\ =>- [(6, r) does not change at t] 

almost sure and vice versa, (6, r) is conditionally constant given events which cause changes 
of (a,q) (and vice versa) so the conditional distribution of £[i]|£[i-] uniquely determined by 
Corollary 12.41 and a formula for (6, r). 
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In the rest of the paper, we shall deal only with the ask side, i.e., with ( A , a, q ), leaving (B, b, q) 
aside due to the symmetry. 

Because density hi depends on all the parameters related to A (i.e., n,p, 6 , t), it is straight¬ 
forward to estimate these parameters by Maximum Likelihood based on h t and a sample from 
(Af w ,e[i],a w ,g[i]). Moreover, if we put 

i = 0 ( 5 ) 

and, quite realistically, assume that k and p depend only on a distance to the ask, i.e., 

n(a,b,p) = k(a- p), p(a, b, p) = p(a - p) (6) 

for some k and p, and that the “inside-the-book” intensities 

k|{l,2,...}, P\ (1,2,...}) (7) 

depend on different parameters than the “spread” intensities 

^|{0 ,-1,-2,...}) p|{0,—1,—2,...}, (8) 

then the estimation of the “ask-side” parameters could be split into the estimation of the 
parameters underlying (JHJ) , based on /j(Aij, e*)^, and that of the parameters underlying (??), 
based on q^ieN- In the present paper, only the latter case of an estimation is discussed. 


3 Generalised Model 

Despite the popularity of the ZI models, there is no doubt that they are too rough and neglect 
many aspects of real-life trading. The most obvious issues in this respect are two: 

(i) the volumes are non-unit in reality, 

(ii) agents, at least sometimes, act strategically rather than in a random way. 

To deal with (i), authors of ZI models usually argue that the sizes of a majority of orders 
are one-lot and, when an order is larger, it may be imagined that several subsequent one-lot 
orders have been issued. However, this simplification may be tolerable only in the case of limit 
orders (if the order book is observed only at certain time instants then larger orders are indeed 
interchangeable with several subsequent unit ones) but it is unacceptable for non-unit market 
orders which could not be imaginatively split without a serious violation of the assumption that 
their arrival intensity is constant. 

To justify the ZI models against (ii), it is usually argued that the strategic reasoning is so 
complex that taking it as random is less evil than constructing wrong models. This argument, 
however, has a drawback, too: no matter how reasonable it may sound, one still has to deal with 
empirical phenomena stemming from (bounded) rationality, such as shifts and rapid insertions 
and cancellations of orders (especially quotes) in response to changes of the book. 

To meet those issues, the setting defined in Section 12.11 might be generalised by assuming 
tbat 
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M volumes of market orders are possibly non-unit, 

51 once a (unit) limit order stops to be a quote due to an in-spread sell limit order, it is, with 

probability 1 — 77 , immediately cancelled or shifted to the position of the new quote, 

52 once a quote jumps out of the spread, there is a non-zero probability that the move was 

caused not by a market order or cancellation but by a shift of the quote. 

Moreover, to give more freedom to modelling of (possibly complex) behaviour of the quotes, it 
is allowed that 

D the distribution of Atj, ej|£[j_i] is arbitrary and, moreover, q and r may jump by more than 
one unit. 


The list of events potentially changing x at ti is newly as follows: 


code 
BMO(z) 
SLO(p, z ) 

CA(a[j_i], *) 
SAL(^, a)] 
SAR(z, a)] 


description 

a buy market order of volume z 

a sell limit order with price %-i] < p < a[i-i] and size z put into the spread 
a cancellation of 2 orders of the ask 

a shift of z orders with price ap_i] to tick ap_i] — a, a > 0 
a shift of z orders with price ap_ 1 ] to tick ap_ 1 ] + a, a > 0 


plus symmetric events concerning b. 

Even though, given this generalisation, (i) of Theorem 12.21 ceases to be true, its Assertions 
(ii) and (iii) would keep holding with r\w v ^ instead of w v ^ if it is additionally assumed that 


• the volumes of market orders are i.i.d. random, independent of all the past and the 
present events on the market 

• the (Bernoulli) variables indicating the shifts are mutually independent, independent of 
the past and the present events on the market H 


What does change, however, is the distribution of quote jumps outside of the spread after the 
ask is depleted, which are now non-zero if and only if 


[ei G {BMO(z), CA(< 7 [i_i]), SAR(z, a)}, z > q[i-i],a > 0], 


(9) 


in which case 


(a[i],?[i]) = (a,q) 44 


Ma,i A Si, M a ,i + 
-1A • 1 ./ A Sj 


Si + q, a <n, 

a — n + 1, q — 0, 


where 

a—l 

= E 4-1 

J =a [i-1] + 1 


3 To see it, note that the proof continues to be valid with P(Bi(g[j_ 1 i,?;) = a) instead of l[a = in (1X1) . 
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is the total number of the orders in the book between the former and the present value of the 
ask and 

i] if e% = BMO(z), 

Si = < — z if e* = SAR(z, a) for some a, 

[o if &i = CA(a [i _ 1 ]). 

is the number of the orders the particular change would like to “eat” from the inside of the 
book. 

As, by (ii) and (iii) of (the modified version of) Theorem 12.21 


M a ,i\£[i-] ~ Bi(z/ a[ ._ 1 ]+ i i j, i]w ai _ 1+ \ : i) o • • • o Bi(r / a _i ) j, rjw a -i t i) o Po 


a— 1 

E 

U-“[t-l] + 1 


l e j,i + L j,i] 


and, by (iii) of the Theorem, is conditionally independent of M a/l given £[,_], we are getting 
Corollary 3.1. On set (0J), the distribution of (a^j, ) I<6 [i— ] is given by density 

~ ( | f s = I EjLo F l M i = j\^[i-M A [i-] = s i + v - j\£[i-]) a [i- 1] <a<n, 

[t 1 \P[Mi < Si|^[i_]] a — n + l,q — 0. 

Remark 2 . As both the value of the ask and its volume are uniquely determined by outside 
set dU]), the distribution of («[*], %;] )!£[*-] is Dirac outside set dHJ). 

Thus, if we keep assuming (J5]) and (J 6 J) , we may estimate parameters underlying /c| {i, 2] ...} and 
y51 {i 2 ,...} together with rj by means of MLE based on gf, as it could be checked in Appendix iBl 
the asymptotic properties of the estimator are not harmed by the generalisation provided that 
the generalised process is ergodic, which may be guaranteed, e.g., by requiring that 


E qt + r t is stochastically dominated by a Markov chain with a zero recurrent state. 

Under this assumption, Proposition 12.11 keeps holding because none of the generalisations, 
except for D, which is treated by E, increase the total number of orders in 5 in comparison 
with the ZI. 

Remark 3. Unlike in the unit-volume Zl model, price impact may be predicted in the GZI 
model: denoting m; the impact of a trade with volume z at time t, we have 


[rrii(z) = m|£[»_i] = £, U = t, , e* = BMO(z)] = P[Ao[i] = = £, U = t, e* = BMO(z)] 

OO 

= 9i(a[i~i] + m, g|f[i_i] =£,U = t, e t = BMO(z)). 

5=0 


4 Empirical Evidence 

4.1 Data 

For our empirical analysis, tick-by-tick trade and quote data was used, i.e., values 

(•£[*] i T*])*gn 
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where is an amount traded at L, defined by 


T[i] 


z ei = SMO(z) 
-2 e* = BMO(z) 
0 otherwise. 


In particular, data of six stocks 

• Exxon Mobile (XOM) 

• Microsoft (MSFT), 

• General Electric (GE) 

• MarketAxess Holdings (MKTX) 

• J2 Global (JCOM) 

• American Realty Investors (ARL) 
from nine US electronic markets 

• NASDAQ OMX BX 

• NSE 

• Chicago 

• NYSE 

• ARCA 

• NASDAQ T 

• CBOE 

• BATS 

• ISE 

from ten months starting from March 2009 were analysed. 

Our choice of the titles was done so as to cover both “large” and the “smaller” stocks. The 
first three ones - XOM, MSFT and GE - belong to the top ten companies by market capital¬ 
ization and usually exhibit large trading volumes. MKTX and JCOM, on the at her hand, are 
usually ranked as “small caps” (having market capitalization around $1 bilion) with moderate 
trading volumes. Company ARL, belonging to the category of “micro caps” (capitalisation 
around $100 million), could be counted among illiquid stocks. 

As to the markets, we did not pre-select which of them to analyse; contrarily, all the markets 
were included on which the examined stocks were traded and for which at least some data were 
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available^- the only exception is NASDAQ ADF, which is a platform for recording trades rather 
than a limit order market (see FINRA (2015))). 

As the trade data (process r) came from a different data source than the quote data (process 
x), it was necessary to algorithmically match individual trades with the LI data changes in 
which our algorithm, designed for that purpose, succeeded in about 70 percent of trades (the 
rest could not be uniquely attributed to any LI change ). This fact, however, did not harm the 
estimation other way than by decreasing the sample sizeJl 


4.2 Estimation 

For the actual estimation, only records originating between 9:40 a.m. and 3:30 p.m. were 
used when the process could be assumed to be near its stationary distribution; the inclusion 
of the ten-minute “warm-up” period following the opening at 9:30 also partially justifies our 
assumption ([5]). 

When estimating within the ZI models, all the jumps of a up@ could serve as a sample 
while, in the case of the GZI models, only the jumps of a matched with trades were included; 
the reason for this restriction is that the jumps of a caused by unpaired trades, shifts and 
cancellations could not be distinguished from each other in the LI data, so the corresponding 
values of s, could not be determined. For each paired trade record, on the other hand, the 
value of Si may be obtained from the volume of the trade and from the (known) volume of the 
former ask. 

For both the ZI and the GZI settings, and for each stock-market pair, the following estima¬ 
tion procedure was run: as its first step, a simple version of the model with 

«(•) = «o, p(«) = po, p 0 > 0, k 0 > 0 (S) 

in its ZI variant), which we later call basic model , was estimated. Subse¬ 
quently, three power tail models with k and p, defined by 

i i < n 

K n _i(i — n)“ K i>n 

were estimated gradually. As a “true” model, the one was chosen in which all the parameters 
came out significant while the likelihood ratio comparing this model with the subsequent one 
came out insignificant. 

Each estimation procedure was run on a sample of at most 5, 000 observations; even though 
there were many more observations available for a majority of the stock-market pairs, a restric¬ 
tion had to be made due to the large time requirement of the estimation, caused by complex 
formulas for conditional densities involved. In case that the optimisation algorithm could not 
find maximum of the ML function within a time limit of 2, 500 sec (approx. 42 min) for none 

4 Our dataset came from Tickdata, Inc. 

5 To be absolutely rigorous, we should assume that the success in matching is stochastically independent of 
the quote process to claim this. 

6 The jumps of a down and the changes of (a;, r) preserving a could be omitted without any loss of information 
because the corresponding these changes are conditionally constant hence not depending on any parameter. 



~( J Pi— 1 i — n (rp \ 

P{V = S x a . ( Tn) 

p n _i(t-n) Qp i > n, 


n = 1, 2, 3, 


( S m ith et al. (20031 ) 
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of the four variants of the model, the procedure was repeated once more with a sample of size 

1 , 000 . 

The prediction power of each model m was evaluated by 

„ , ES M i4i- E (4]i^-])i 

* m 1 N+M I + ~i_ | 

Ej=w I a [i\ a + N\ 

where N is a sample size, M = 1V/10 is the number of out-of-sample observations, af ,..., cl^+n 
are the magnitudes of the jumps of a up, is their average and where the conditional expecta¬ 
tion in the numerator was computed by means of the estimated parameters. By its construction, 
P m can be understood as a percentage improvement in comparison with a naive prediction of 
the out-of-spread jumps by their mean. In the GZI setting, P m may also be seen as a measure 
of accuracy of the market impact prediction. 

For automation of the whole procedure including matching of trades, and the model se¬ 
lection, a C++ package has been developed by the author, using the NLOPT library, namely 
LBFGS and BOBYQUA optimisation algorithms, to compute maxima of the ML functions Jj 
A source code of the package is available at https://github.com/cyberklezmer/fepp under 
branch qfl5. 


4.3 Results 


Results of the estimation procedure of the ZI modelsd for the “big-cap” stocks XOM, MSFT 
and GE are summarised in Tabled] It may be seen that 21 out of all the 27 estimations were 
successful in the sense that a significant model was found, which exhibits a positive improvement 
(measured by P m ). Out of the “failed” cases, four times a significant model was found giving a 
negative -P, r jj, once the time limit was reached even with a reduced sample, and once no data 
was available. 

Out of all the 21 successful instances, the basic model (S') was chosen nine times, the simplest 
tail model (T]) six times and model (T 2 ) six times; model (T 3 ) never came out significant within 
the time limit. Higher P m ’s values in comparison with a basic model were reached in only seven 
out of the 12 cases when a tail model was selected, while the basic model came out better in 
two cases; in the three remaining cases, the P m was identical up to two decimal digits. Even 
though, on average, the tail model performed better, the difference was not found statistically 
significant according to a Wilcoxon test. Detailed results of this estimation can be found in 
Appendix 1C. 11 

Summary results for the GZI model for big-caps may be found in Table [2j Here, the 
estimation was successful in 21 out of 25 cases with at least some data, three times the time¬ 
out was reached; P m came out negative in one case. 

In order to answer the question whether it is worth using the more complicated GZI model 
rather than the ZI one, the estimation procedure was run for ZI models once more using the 


' As a default choice LBFGS has been used, the BOBYQUA was employed only in instances in which it could 
not find a maximum within the time limit). For more on these algorithms and the NLOPT library, see .Tohnsonl 

{ml). 

8 In two of these instances, all the jump increments were exactly one; hence their average predicted the jumps 
with zero error, which lead to minus infinity for P m . 
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same samples as those used for the GZI ones (i.e., only the changes with trades matched); the 
corresponding results are displayed in the second column below each stock. Out of the 20 cases 
when both estimations were successful, eight times the GZI model reached a higher P m while 
three times P m was higher for the Z1 model; nine times P m was the same up to two decimal 
digits. Again, it could not be statistically proved that the GZI model performs significantly 
better. 

Detailed results of the GZI estimation may be found in Appendix IG.2I (the GZI models) 
and Appendix IG.3I (their Z1 counterparts). 

Tables [3] and [4] summarise estimates of tail exponents of k and p. It may be seen that a 
dispersion of the values of a K is great; h owever , their average —2.3(0.41) is not far from present 


empirical evidence (see, e.g., Chakraborti et al. ( 2011a ). III.C)@The results for tail exponents 


of p, which have never been statistically estimated yet to the best of the author’s knowledge, 
are similar to the case of a K , with average value —2.48(0.42). 

Results of the estimation of the small caps MKTX, JCOM and the micro-cap ARL may be 
found in Tableland Table [6j For the small caps, the ZI model was successful for 10 of 13 pairs 
with sample size at least IV > 20, once a negative P m was reached, two times the optimization 
failed@ For ARL, the procedure was successful only once (with a disappointing prediction 
power) with all of the failures due to failed optimizations. A closer look to the results suggests 
two possible causes of the failures: small sample sizes and/or large average jumps of a, causing 
large evaluation times of the densities (see Appendix O for detailed values). The situation is not 
much better in the case of GZI models either: small caps are successful in 7 out of 9 cases, ARL 
is only once half-way successful out of three cases. Results of a comparison of more complex 
ZI models with their basic variant are similar to those of the big-caps: out of 8 cases when the 
comparison is possible, more complex models won three times, four times the basic model was 
more successful, once the results were the same. Similarly, from 8 comparisons of GZI and ZI 
variants of the model, the GZI variant was more successful four times, the ZI one three times, 
once the results were equal. 


5 Conclusions 

A setting covering many of the existing zero intelligence (ZI) models and its generalisation 
allowing for non-unit market orders and shifts of quotes (GZI) were introduced. Several variants 
of both the ZI and GZI settings, differing in their complexity, were tested on trade and quote 
data for 54 real-life stock-market pairs. It was found that, especially for liquid stocks, both 
the ZI and GZI models came out significant with a substantial prediction power; however, 
their more complicated variants did not produce significantly better predictions in comparison 
with a simple model with constant intensities. Finally, as the rewults of the ZI and the GZI 
variants are comparable, our suggestion is to use the GZI models which are more realistic and 
are capable of predictions of market impact. 

9 It should be noted, however, that only the first several values of k and p played a role in the estimation, 
especially when the average jump out of the spread was small, so the results say very little about the actual 
power law behaviour. 

10 In particular, the solver stopped the optimization without changing the initial parameters reporting that 
its tolerance criterion was reached which suggest that the MLE function is too flat to be optimized. 


15 











XOM 

MSFT 

GE 

NASDAQ OB 

Ti(< 0) 

77(0.70,0.70) 

5(0.48) 

NSE 

77(0.91, 0.90) 

77(0.89, 0.87) 

Ti(0.76, 0.76) 

Chicago 

77(0.47, n/a)~ 

5(0.93) 

T 1 !(0.76, 0.71) 

NYSE 

5(0.33) 

X 

5(0.53) 

ARCA 

Ti (0.92, 0.91) 

5(—oo)* 

t 

NASDAQ T 

7i (0.53, 0.53) 

5(1.00)* 

5(—oo)* 

CBOE 

S(< 0) 

77(0.35,0.36) 

5(0.67) 

BATS 

77(0.86, 0.83) 

5(1.00) 

5(0.00)* 

ISE 

77(0.66, 0.67) 

5(0.66) 

77(0.94,0.93)* 


Table 1: Results for ZI. Notation: m(p,p 0 ) - model m was selected with P m = p and with 
Ps = Po, minus sign - sample less than 5000 was available, * - a reduced sample used, t - time 
limit reached when without finding a result, x - less than 20 observations were available. 



XOM 

MSFT 


GE 


GZI 

ZI 

GZI 

ZI 

GZI 

ZI 

NASDAQ OB 

77(0.19) 

77(0.17) 

77(0.68) 

77(0.66) 

Ti(0.74) 

5(0.74) 

NSE 

77(0.70)* 

77(0.69)* 

77(0.38) 

77(0.29) 

5(0.33) 

t 

Chicago 

77(0.35)- 

77(0.34)- 

77(0.49)- 

77(0.51) 

^(0.79) 

77(0.79) 

NYSE 

5(0.17) 

5(0.13) 

X 

X 

5(1.00)* 

5(1.00)* 

ARCA 

77(0.84) 

T 2 (0.84) 

t 

5(0.00)* 

t 

t 

NASDAQ T 

5(0.74) 

Ti(0.74) 

X 

X 

t 

t 

CBOE 

77(<0)" 

77 (< 0 )- 

T 2 (0.72)- 

5(0.68)- 

77(0.23) 

77(0.23) 

BATS 

77(0.31) 

Ti(0.60) 

5(0.00)* 

5(0.00)* 

5(0.00)* 

5(0.00)* 

ISE 

77(0.44)* 

5(0.42)* 

T^O.96)* 

T^O.96)* 

5(0.36)* 

5(0.42)* 


Table 2: Results for GZI. Notation: m{p) model m was used with P m = p, minus sign - sample 
with less than 5, 000 available, * - a reduced sample used, t - time limit reached without hireling 
a result, x - less than 20 observations were available 
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XOM 

ZI 

GZI 

MSFT 

ZI 

GZI 

GE 

ZI 

GZI 

NASDAQ OB 

-5.49 

-2.60 

-0.38 

-0.48 

n/a 

-4.15 

NSE 

-1.11 

-2.83 

-1.30 

-2.64 

-0.12 

n/a 

Chicago 

-2.31 

-2.28 

n/a 

-0.75 

-9.72 

-1.47 

NYSE 

n/a 

n/a 

X 

X 

n/a 

n/a 

ARCA 

-2.18 

-1.75 

n/a 

t 

t 

t 

NASDAQ T 

-0.84 

n/a 

n/a 

X 

n/a 

t 

CBOE 

n/a 

-1.70 

-0.97 

-0.34 

n/a 

-2.12 

BATS 

-1.15 

-4.84 

n/a 

n/a 

n/a 

n/a 

ISE 

-1.18 

-2.76 

n/a 

-6.76 

-0.26 

n/a 


Tabic 3: Estimated values of a K . Notation: n/a - only model 5 came out significant, t - time 
limit reached, x - no data available. 



XOM 

ZI 

GZI 

MSFT 

ZI 

GZI 

GE 

ZI 

GZI 

NASDAQ OB 

-5.37 

-2.60 

-0.26 

-0.60 

n/a 

-4.22 

NSE 

-1.84 

-2.89 

-0.89 

-2.31 

-0.61 

n/a 

Chicago 

-2.13 

-3.22 

n/a 

-1.94 

-8.28 

-1.77 

NYSE 

n/a 

n/a 

X 

X 

n/a 

n/a 

ARCA 

-1.77 

-0.77 

n/a 

t 

t 

t 

NASDAQ T 

-0.67 

n/a 

n/a 

X 

n/a 

t 

CBOE 

n/a 

-2.21 

-0.03 

-9.24 

n/a 

-3.39 

BATS 

-0.38 

-4.59 

n/a 

n/a 

n/a 

n/a 

ISE 

-0.94 

-2.43 

n/a 

-3.56 

-0.66 

n/a 


Table 4: Estimated values of a p . Notation: n/a - only model 5 came out significant, t - time 
limit reached, x - no data available. 



MKTX 

JCOM 

ARL 

NASDAQ OB 

r 

7i (0.06, 0.06) 

X 

NSE 

X 

r 

X 

Chicago 

X 

X 

X 

NYSE 

X 

X 

/ 

ARCA 

T 2 (0.40, 0.45) 

T 2 (0.18,0.17) 

/* 

NASDAQ T 

T 2 (0.44, 0.46) 

T 2 (0.38, 0.37) 

5(0.08)' 

CBOE 

5(0.92)* 

t 2 (< 0, < 0) 

r 

BATS 

T 2 (0.13, 0.35) 

T 2 (0.60, 0.52) 

f~ 

ISE 

5(0.31)* 

T 2 (0.39, 0.49) 

f~ 


Table 5: Results for ZI of the small caps. Notation: / - optimization failed, for the rest, see 
Table [0 
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MKTX 

GZI ZI 

JCOM 

GZI ZI 

GZI 

ARL 

ZI 

NASDAQ OB 

X 

X 

T 2 (0.78)- 

Ti(0.79) - 

X 

X 

NSE 

X 

X 

X 

X 

X 

X 

Chicago 

X 

X 

X 

X 

X 

X 

NYSE 

X 

X 

X 

X 

5(< 0)- 

- t 2 {< 0)- 

ARCA 

T 2 (0.31)- 

T 2 (0.25)- 

5(0.08) 

T 2 (0.04) 

f~ 


NASDAQ T 

X 

X 

X 

X 

5(< 0)- 

- Xi(0.49) _ 

CBOE 

5(1.00)- 

5(1.00)“ 

71(0.29)- 

Xi(0.25)- 

X 

X 

BATS 

5(0.41)“ 

5(0.39)- 

T 2 (0.65) 

5(0.66) 

X 

X 

ISE 

s(< or 

s(<oy 

5(< 0)- 

5(< 0)- 

X 

X 


Tabic 6: Results for GZI of the small caps. Notation: f - optimization failed, for the rest, see 

Tabic m 
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A 


Proof of Theorem 


12/2 


Before starting the proof, recall that that a Markov chain X on M + is called immigration and 
death process with immigration intensity k and death intensity p and with initial population 
A (we abbreviate this by Z(n,p,A)) if its transition matrix A = (Ajj) has zero components 
except of \jj+i = k, Xjj-i = jp , j > 0 and if X 0 = A. Notice also that number of customers 
in an M/M/oo in queuing theo ry model follows an I process. It is well known (see, e.g., 
iMandies and Zuraniewsld ( 201lh . Section 2) that, for any positive t, 


k , 


l(K,p,A) t ~ Po ( -[1-exp {-pt}} ) o Bi(A, exp{—ps}). 


( 10 ) 


Further, let us introduce the following (re)formulation of process S restricted to time interval 

(0,ti], which will be used repeatedly in the subsequent proof: Let u 1 , y 1 , u 2 , v 2 . u n , v n be 

independent uniform variables, independent of S 0 . By Kallenberg (2002), Theorem 6.10, there 
exist mappings U 1 , Id 1 ,..., f/ n , V n from Ng x [0,1] to the space of stochastic processes on M + 
such that, for each 7T, U n (A 0 , u 7r )|S 0 ~ I(^ |0 ,^,o,^o) and V *(B 0 , v n )\-o ~ cr nfi , B%). 

By Kallenbera (2002), Proposition 6.13., U 1 , V 1 , U 2 , V 2 ,..., U n , V n are mutually conditionally 
independent given S 0 , so we may assume, without a change of distributions, 


At = (Ul, • • •, un, B = (V t \ ..., V t n on [t < h}. 


( 11 ) 


(to see it, check definitions (JT|) and (J2J)). 

Using this reformulation and noting that both e 1? At\ are functions of ( U 77 , U 7r ), r < ao , we see 
that 

(t uei )AL So (Un,> ai (12) 

where AALcB denotes conditional indepdence of A and B given C. 

Because, by the definition of the ask, = • • • = 1 = 0 and A °°_j = qo, i.e, 

..., A^n are conditionally constant, we have, for any si, s 2 ,..., s n , 


IP[^4[i_] — si, — s 2 ,.. •, — s n |ei, U, s 0 ] 

= l[si = s 2 = ■ • • = s ao _i = 0, s ao = qo] • P [U^ 0 - 1 = s a o+i, • • •, = s n |ei, U, So] 

n 

= l[si = s 2 = • • • = s ao _i = 0, s ao = qo}- JJ p [U/ = s/ (13) 

7T=ao+l 


— (ni a oTl /i/Cto+2 


where the last follows from Hoffmann-Jor g enson (1994). (6.8.14) (with R = (u'" u | i , u“ 
and S = (ei,ti,S 0 )) and from the fact that the jumps of f/’s do not coincide with probability 
one (which guarantees that U£ — U£ , n > a 0 , almost sure). 

As a first step of proving the assertions of the Theorem, let us deal with one-step-ahead 
distribution P[(A ti,e-i,A[i-\) £ • |S[ i _ 1 ]]: 


u 


Proposition A.l. (%) For any r and e, 


P[A U > r , e-i = e|H[j_i]] = exp{-7j_iT}7Tj(e). 
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(ii) For any p, 


|Ci) Afj, 1 —'[i—1] 


(S 0 

l ^i[i~ i] 

[po (0p,j-i(l — ^p,j-i)) ° Bi , 5p,i-i^J 


on [1 < p < a[j_i]] 
on \p = a[j_ij] 
on [a[j_i] <p <n\. 


(in) For any p, 


A«] Aij+l, *£[j]j At j, 1 —'[i—1] 

J<5a h on [1 < p < max(a [j ],a [i _i])] 

\ Po (0 Pl[ i_i](l - 5 p ,i_i)) o Bi (Ajj-ip <Jp,i_i) on [max(a M , a[i_i]) < p < n]. 


Proof. Thanks to the homogeneity of the process, we may assume i — 1. 

(i) It follows from textbook knowledge that Ati, being the first jump time of Markov chain 
(V bo ,U bo+1 ,V ao+1 ,...,U ao ), is exponential with rate 70 , while e*, coding a type of the chain’s 
first jump, is (conditionally w.r.t. S 0 ) independent of A t\ with probabilities of particular events 
being equal to the rates of the events’ intensities to 70 , which is formally expressed by (i). 

(ii) The formula follows from (fT3]) and ffTUl) . 

(iii) Let 1 < j,k <n and put Sj ik = [np] = j, o,q = k\. Clearly, 


£[i] = X[ij on S jjk , (14) 

where 

%] = (j, A J {1 j, 6 , B b {1] ), b = max{ 7 r < j : B( {] > 0}. 

Further, it follows from the dynamics of the process and from the almost sure exclusivity of 
jumps of 5 that 


Ah = ^ on S jtk , A = (r jtk , A$*+\ Af H ), 


(0,..., 0, gpj) 


j > k 


r j,k = 


j* n 5 ) 

(O^^O, q m , , go) j < k. 

jx (k-j- l)x 


Now, because (Ai 1; xjij) are functions of U 1 , V 1 ,... , S 0 , we have (Ati, xq]) G cr(t/ 1 , V 1 ,..., W, , 

so we may use the Law of Iterated Expectation to get 


P[i G .|%], Aii, S 0 ] = E(P[i G *\U\ V\ ..., U\ V j , S 0 p[i], A t u S 0 ) 

= E(P [A G •|H 0 ]|x[i ] , Afi,S 0 ) = P [A G «|So] = p jt *,(•; B, E 0 , r jtk ) on S jik , 

n 

Pj, k {s\, ■ ■ ■ ,s n ;i,S,r) = l[ri = Si,. ..,r jvk = s jyk ] F[U n (E) t = s*]. (16) 

n=j\/k-\-l 
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(at the third we have used Proposition 6.6 of Kaljenberg (120021 ). the fourth one is due to 
measurability of the inner term w.r.t. the condition of the expectation, the fifth one is due to 

(USD). 

Now, by Local Property (IKallenbergl ( 20021 ). Lemma 6.2) applied to ffTLl) . (1I5|) and (fT6l) , we 
are getting that 


P[A[ 1 ] G •|X[ 1 ], Afi, £1 0 ] — Po [ 1 ] ,ao(*; h, “0, (^[i], jtfq, . . . , Aj 1 , 11 °)) 


(17) 


on each Sj^ and, as the Sj t k s cover all the underlying probability space, the relation holds 
universally. 

Finally, by the strong Markov property of S (guaranteed, e.g., by iKallenbergj (120021 ). The¬ 
orem 12.14.), we get 

P[A t 2 > r,e 2 = e|S[i],X[i], Ati,E 0 ] = P[At 2 > r, e 2 = e|E[i]] 


which is, by Proposition 6.6. of Ka llenber g (2002j), equivalent to 

A t 2 , e 2 -LL X[1]) Ati,s 0 “[i] 

implying, by the same Proposition with switched variables, 


P|A[i] G *\ e 2, At 2 , X[i], Ati, ^o] — P|A[i] £ •|^[i] ) At l5 ^o] 

which, together with (UTl) . proves (iii). □ 

Proof of the Theorem 
Ad (i) 

By the Law of Iterated Expectation, the strong Markov property of 5, and Proposition lA.il (i). 
it follows 


P[A U > T,ei = e|£[i_i]] = E[P[Atj > r, e* = e|E[j_i], 

= E[P[Atj > T 1 e i = e|S[i_i]]|^[i_i]] = E(exp{-7 i _ir}7r i (e)|^_i ] ) = exp{-7j_ir}7Ti(e) (18) 

(we could get rid of the outer conditional expectation at the last equality because its integrand 
is measurable with respect to its condition). 

Ad (ii) 

Before deriving the distribution of A[j_], let us formulate an auxiliary result: 

Lemma A.2. Let A, B be random variables, let T be as a-algebra and let S G T. If B\T ~ 
Po(e) o ~Bi{y,zu) on S and A\B, T ~ Po(£) o Bi(.B,<5) on S for some T-measurable variables 
e, is, w, C, 5 then A\T ~ Po(£ + 5e) o Bi(^, Szu) on S. 
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Proof. Let A 0 , B 0 , C\, _D 1; Ei, C 2 , D 2 , E 2 ... be variables mutually conditionally independent 
given T , such that, on S, AqIJ 7 ~ Po(£), B 0 \E ~ Po(e), and, for each i G N, CfT ~ 
Bernoulli(5), DfJ 7 ~ Bernoulli^), ~ Bernoulli (m) (outside S, all the variables could 

be e.g., zero). Assume, without a change of the examined distributions, that B = Bq + E\ + 
E 2 -\ - \-E v and A = A 0 + A 1 + A 2 where A 1 = Y,f=i Ci, A 2 = XX=i DjEj. Clearly , A 2 is (condi¬ 

tionally) Binomial with parameters v and 5m. Further, by Dalev and Vere-Jonesl ( 2003 1. 2.3.4, 
Ai\E ~ Po(he). As Aq, A x , A 2 are mutually conditionally independent given T (which is because 
each of them depends on independent variables) we get that A = Po(() o Po(<5e) o Bi(z /,5m). 
The Lemma now follows from the fact that a convolution of two Poisson variables is again 
Poisson. □ 

Let 1 < p < n. First we prove, by induction, that, for any k, 

A[ fe ]|£[( fc+ i)-] ~ Po(e Pi fc + t Pjfe ) o Bi(u Pjk , m P)k ), on [a [fc ] < p < n\. (19) 

First, we show that (fl9ll holds for k — 0: indeed, by Proposition IA.1I (i) and Kalle liberal 


(120021 1. Proposition 6.6., we get that (ti, ei)_LL a;o So (recall that £[i_] = xo), which, by the 
Proposition 6.6. and by ([3]) , gives (Toll for k — 0. 

Now, assume (fT9l) to hold for k — 1. Similarly as in (IT8lh we get, by Proposition I A. 11 (iii), 
that 


P [^[fc] G # l^ffc-i]^[(fc+i)-]] 

= E(P[Af fe] G •I4- I ]]l4- I p W)-]) = E ( P ^fA 


Ik] 


= P 


- , [fe—i] j ^k, Atfc] | A[^_i], £[&_]) 


Po (lb. 


giving, by Lemma lA.21 applied to the induction hypothesis and (|20jl . validity of (1T9T1 on [a^-i] < 
p < n\. As, on [a[j_p = p\, AC = q^} is conditionally constant, i.e., 

p [^ffc] = a lC[(fc+i)-]] = 1[« = ml = p [ Bi (^f fc p 1) ° p o(0) = «] 

(fl9|l is proved on [a[i—i] = p\ because k k (p) = k on the set so e P)k = i P , k = 0, m Pjk = 1 and 
u Pjk = AC on the set. As (fl9|l on [ap_i] > p\ follows similarly (here, AC = 0), (fl9|l is proved 
for k. 

Having proved ([H?]) for any k, we may proceed to (ii) of the Theorem: Let 1 < p < n. 
Similarly as above we get, using Proposition IA.1I (ii). 


PP?H 6 •I4- 1 iX[h] = EtPPfi. 


-l]) e.j]|A|’_ 1 j, ^[j_]) 


= P 


Po (0 p ,i-i(l - 6 p ,i-i)) o Bi (A p {i l] ,5 p ^A G • on [a[j_i] < p < n] 


which implies (ii) on [a.p_i] < p < n] by Lemma [A.21 Validity of (ii) on set [p < a[*_y] on which 
Aj)_j is conditionally constant, follows similarly as in the proof of (fT9]h 
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Ad (iii) 

Let i G N, 1 < p < n, denote A~" p = (A^^p, and assume that are conditionally 

independent given %_], i.e., that 

P[4-,l = r, A^_ n = RK tt . } = «] = P[4_„ = r|« IH = = 5] (21) 

holds for any r, R , £ and for k = i. 

Let / C Mq and put 

q(s, e, t, S) = Pe /|ei = e, At ± = r, S 0 = S], 

From (fT3]h we have that 

q(I,e,T, So) = JI q v (I v ,e,T,A%) 

{l<TT<n,I n ^ 0 } 

for some g 1; ..., q n . So, by the Law of Iterated Expectation, strong Markov property and the 
homogeneity, applied gradually, 

P[A[j_] e I\£[i~], = E(P[A[j_] e /|^[i-],5[ i _ 1 ]]|^[ i _], A[j_!]) 

= E(P[A[i_] e I\ei, Ati, S[j_i]]|£[j_], A[j_i]) 

= E(q(I,e i ,At i: E [i _ 1] )\^ [i _ ] ,A [i - 1] )= q n (I„, e i; AR, Aj^) (22) 

{l^n^No} 

(we could get rid of the expectation on the RHS because its integrand is measurable). 

Let 1 < p < n, let J C N 0 and K C Nq _ 1 . By evaluating (|22|) with / = N 0 x---xJx---xN 0 
and with / = K 0 x • ■ • x No x • • • x Ji p , we get 

p [4-l e e VijI 

= pyy, e e x|« [H , Vij] 

= PI4-1 e e 

(we could write the last “=” because the right-hand side of (l22|) does not depend on no A n 
with I n = No). Thus, by the Complete Probability Theorem, 

= E E P Nh e J ' T-] e ^1^5-11 = A-H = R ’ %-l = a 

r _R 

x p N,_ 1| =''.2l | l p _ 1| = flk|H=«] 

E P Nh e J ITi-H = r ^[H = f] p l4-i] = r l«IH = a 

r 

* (EW-, e K AA = «.«:<-! = = «l%-i = £lj 

= PP5-, e J|«[H = «]P[y', e = {) (23) 
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proving conditional independence of Aj*_j of Ajl^ given £[j_] for any p, which implies mutual 
conditional independence of A^, A|_j, • • •, Aj^_] given £[*_] (see Kallenbergl ( 2002 1. p. 109 
below). 

ft remains to prove (l2Tjl for k — i + 1. To this end, put Sj,k,r,s = [ap] = j,a[i- i] = k,q\{\ = 
r, q\i-\] = s] for any j, k, and observe that, on each Sj^, r ,s, Aj^,..., A J A k are conditionally con¬ 
stant while Ajq = Ajj_j, for n > j \Jk. As A^,..., Aj*_j are mutually conditionally independent 
given £[j] by (123|) and as (conditional) constants are trivially (conditionally) independent of any 
random element, we get that Ah, ..., AT are mutually conditionally independent given on 
Sj,k,r,s by Local Property ( Kallenbergl ( 120021 ). Lemma 6 . 2 ), and, as S’s cover all the underlying 
space, the conditional independence holds universally. Relation (l2Tjl now follows from (i) of 
Proposition IA. II (the same way as in proof of (iii) of the Proposition). 

Finally, as (I2T|) holds for k = 1 by (J3J), we have just proved (iii) by induction. 


B MLE Estimation 


In the present Section, asymptotic properties of the MLE estimator are proved both for the ZI 
and the GZI model. 

Before starting, let us formulate an auxiliary result: 


Lemma B.l. If X is a continuous time stationary ergodic Markov chain in countable space 
S then Y = (Ar*, X r .) ieN , where t 1 , t 2 ,... are the jump times of X, is a stationary ergodic 
stochastic process. 


Proof. Denote A (a;) the intensity of th e first j ump of X given that A" 0 = x. From th e stro ng 
Marko v property (Theorem 12.14 of Kallenbergl ( 2002 )). from Lemma 12.16 of Kallenberg 
(120021 ) and from the scalability of exponential distribution we have that U \, LL,..., where 
Ui = ATj / A(A ri _ 1 ), i G N, is a sequence of i.i.d. (unit exponentially) distributed variables, 

independent of An. X. 


an , 1 


„_, As U n - being an i.i.d. sequence - is a strong mixing and X Tri is a 

strong mixing by Bradlevl ( 2005 1. process Z n := (X Tn , U n ) is a strong mixing (note that, for any 


A = A 1 x A 2 , B = B 1 x B 2 , A 1 , B 1 e S z , A 2 ,B 2 e K z , it holds that P (Z e AnT n B) = P {Z 1 e 
(A 1 n TnB 1 ))P(Z 2 e (A 2 n T n B 2 )) p (z 1 e A 1 ) P(z x e T„R 1 )P(Z 2 e A 2 )P(Z 2 e T n B 2 ) = 
P(Z E A)P(Z E B) where T n is a shift operator; the case of general A, B follows from their 
approximation by rectangles) we get the ergodicity of Z by the well known fact that strong 
mixing implies ergodicity. Finally, as Y n is a function of Z n , the Lemma is proved. □ 


It follows from the Lemma that, once is understood as a process with t E (—oo, oo) then, 
thanks to its ergodicity (see Proposition 12. Ill process £[»_], i E Z, is also ergodic. Hence by 
Birkhoff-Khinchin theorem ( Cornfeld et all (19821), Appendix 3), for any measurable function 
/, 


1 




(24) 


Z— 1 


in probability given that the expression written on the RHS exists. 
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Now, let 9 0 G be a vector of true parameters. Denote 


B n — E [L n (9 0 )] , 


N 

L n (8) = 

2—1 


{ 9. 2 («) 


1 

m=l,...,k,n=l,...,k 


( B n being called Fischer information matrix) where g (m ' ) stands for g. Substituting map¬ 
pings assigning ->■ g { and £[*_] -» V& for / in ([24]), we get, by ([231), tliat 

±L N = M(9 0 ) (25) 

for some non-random matrix M which, being a limit of positive semidehnite matrices, is also 
positive semidehnite. By taking expectations on both sides of (125|) . we further get 


Now put 


and observe that 



M(6 0 ). 


N 

h,m = E v »i°g(*(®)) 

2 — 1 


In{0) = K n {9) - L n {6), 


K n = 


E 


i 

W) 


vhilfi) 


(26) 


As J gi(6) = 1 by the definition of density, any of the integral’s first or second derivatives has 
to be zero, and, in particular, V 2 f g(9) = ( J \7 2 g{9) = 0 (we may interchange the integral and 
derivative in our discrete case) so 

EAA(^o) = j V 2 ^o)<^o) = j V 2 ^(0 O ) = 0 

implying 

~j^^n(9o) M(6 0 ). (27) 

If we now, as usual, assume M(9 0 ) to be regular, then, by (1261) and ([ 27[) . 

-N~ 1/2 B~ 1/2 I n (9 0 ) = [(N- 1 B n )-\-N- 1 I n (9 0 ))] 1 / 2 [N-\-I n (9 0 ))} 1 / 2 M(0 O )“ 1/2 . (28) 

Let H6 1 — 6*011 = S for some 5 and, for any matrix J, denote F)v(J) = — (9 — 9 0 ) t [N~ 1 /' 2 B n 1 ^ 2 J](9 — 
9a). By f[28l) and basic linear algebra, 

lim F N (I N (9 0 )) = (9 - 9a) T M{9a)- 1,2 {9 - 9 0 ) > 5 2 A min (29) 


where A min is the smallest eigenvalue of M(9q) 1 ^ 2 . If, in addition, the parameter space is open 
and both k and p are twice continuously differentiable with respect to all the parameters (which 
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is true for all the versions of the model we use), then -^1^(6) is asymptotically Lipschitz on 
{6> : ||6> — 6» 0 || < so 


lim\F N {I N {0) - I N {6 0 ))\ 

< lim ||0 - 0o||||^“ 1/2 ^ 1/2 ||||/ iV (0) - 7^(00)11110 - 0o|| < K s \\M(6 0 )- 1/2 \\5 3 (30) 

where || • || is a suitable norm and Kg is a Lipschitz constant, further implying, by using the 
fact that |a + 6| > |a| — \b\ together with (1291) and ([29)1 . we get 

lim F N (I N (6)) > limF JV (7 JV (0 o )) - lim | F N (I N (d) - J^(0 O ))| > 8 2 g(8). 


where g(8) = (A min — Kg\\M(6 0 )~ 1 ^ 2 \\8). From the differentiability it follows that Kg —y K as 
8 —> 0 for some K so there has to exist A and q n > 0 such that g(8) > go for all 8 < A. Weak 
consistency now follows from (2.3) of Crowder ( 19T6h with c n = rP for suitable (3 < 1. 

For the asymptotic normality it suffices, by Crowder ( 19761 ) . (4.13) and the considerations 
explained below, that the absolute fc-th moments of the observations are bounded for a certain 
k > 2. This, however, can be easily achieved by determining a large enough constant L and 
excluding from the sample any observation (rc^pT]*]) with E(||(x[j],T[j])|| fc |£[;_]) > L. 


C Detailed Results 

Notation 

N - sample size, Aa + - average jump of a up, q - average market order volume. The number 
of stars stands for signification on levels 0.05, 0.01 and 0.001, respectively. 


C.l ZI Model 


XOM on NASDAQ OB 

XOM on NSE 

XOM on Chicago 

1 V = 5000, a+ = 1.308000, q = 1.000000 

N = 5000, o+ = 48.6565 

N = 119, a+ = 6.9466 

P s = -0.150 Pm = -0.147 

P s = 0.897 P m = 0.910 

P m = 0.466 

k 0 = 568.7488(140.9081)*** 
ex* = -5.4933(0.3325)*** 

PO = 510.4026(126.1866)*** 
a p = -5.3680(0.3256)*** 

k 0 = 0.5002(0.0227)*** 

= 0.0029(0.0007)*** 
ex* = -1.1057(0.0741)*** 

PO = 6.8241(0.1139)*** 

PI = 0.4923(0.0244)*** 
ol p = -1.8369(0.0464)*** 

k 0 = 0.0187(0.0018)*** 

«l = 0.0800(0.0093)*** 
c* K = -2.3140(0.0786)*** 

P0 = 0.0128(0.0013)*** 

Pl = 0.1151(0.0131)*** 
otp = -2.1322(0.0710)*** 

XOM on NYSE 

XOM on ARCA 

XOM on NASDAQ T 

N = 5000, a+ = 1.6307 

N = 5000, a+ = 1.1156 

N = 5000, a+ = 1.0544 

P s = 0.332 

P s = 0.915 Pm = 0.917 

P s = 0.532 Pm = 0.532 

k 0 = 4.6981(0.0163)*** 

PO = 0.9183(0.0029)*** 

k 0 = 50.8499(0.4649)*** 
ot K = -2.1805(0.0213)*** 

PO = 8.1381(0.0680)*** 
ol p = -1.7708(0.0192)*** 

k 0 = 339.0192(6.6030)*** 
o: K = -0.8425(0.1939)*** 

P0 = 52.1695(0.9973)*** 
otp = -0.6747(0.1931)*** 
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XOM on CBOE 

XOM on BATS 

XOM on ISE 

N = 5000, a + = 64.8869 

N = 5000, a + = 1.0255 

N = 5000, a+ = 8.1611 

P s = -0.019 

P s = 0.827 P m = 0.859 

P s = 0.674 P m = 0.656 

KQ = 0.0000(0.0000)*** 

P0 = 0.0000(0.0000)*** 

k 0 = 649.2523(12.7981)*** 

= 7.6678(0.7012)*** 
at K = -1.1529(0.1952)*** 

P0 = 97.1699(1.8800)*** 

Pl = 1.2875(0.1017)*** 
oc p = -0.3755(0.2070)* 

k 0 = 3.6649(0.0790)*** 

«l = 0.1846(0.0062)*** 
ol k = -1.1812(0.0154)*** 

P0 = 8.7439(0.1230)*** 

Pl = 0.6432(0.0117)*** 
ocp = -0.9377(0.0122)*** 

MSFT on NASDAQ OB 

MSFT on NSE 

MSFT on Chicago 

N = 5000, a+ = 1.0820 

N = 5000, a+ = 1.3504 

N = 5000, a+ = 3.6887 

P s = 0.703 P m = 0.700 

P s = 0.866 Pm = 0.889 

P s = 0.929 

k 0 = 4.7416(0.0174)*** 
o: K = -0.3883(0.0325)*** 

P0 = 0.1427(0.0005)*** 
oc p = -0.2643(0.0319)*** 

k 0 = 0.0088(0.0004)*** 

«! = 0.0129(0.0020)*** 
o: K = -1.3047(0.0966)*** 

P0 = 0.0094(0.0002)*** 

Pl = 0.0050(0.0004)*** 

OL p = -0.8890(0.0653)*** 

k 0 = 0.0020(0.0000)*** 

P0 = 0.0007(0.0000)*** 

MSFT on NYSE 

MSFT on ARCA 

MSFT on NASDAQ T 

N = 0, a+ = 0.0000 

N = 1000*, a+ = 1.000000 

N = 1000*, a+ = 1.000909 

No data 

P s = -inf 

P s = 1.000 


k 0 = 1.2963(0.0034)*** 

PO = 0.0105(0.0000)*** 

KQ = 1.1965(0.0028)*** 

P0 = 0.0064(0.0000)*** 

MSFT on CBOE 

MSFT on BATS 

MSFT on ISE 

N = 5000, a + = 19.3225 

N = 5000, a+ = 1.0002 

N = 5000, a + = 1.3325 

P s = 0.359 P m = 0.346 

P s = 1.000 

P s = 0.656 

k 0 = 0.2417(0.0026)*** 

/ti = 0.0606(0.0006)*** 

cx K = -0.9796(0.0050)*** 

P0 = 0.0825(0.0009)*** 

Pl = 0.0168(0.0002)*** 
oc p = -0.0336(0.0055)*** 

K 0 = 1.3741(0.0010)*** 

P0 = 0.0245(0.0000)*** 

k 0 = 1.9125(0.0027)*** 

P0 = 0.0551(0.0001)*** 

GE on NASDAQ OB 

GE on NSE 

GE on Chicago 

N = 5000, a + = 1.0411 

N = 5000, a+ = 1.0415 

N = 5000, a + = 2.8520 

P s = 0.481 

P s = 0.763 P m = 0.757 

P s = 0.714 Pm = 0.758 

k 0 = 1.3310(0.0064)*** 

P0 = 0.0872(0.0003)*** 

k 0 = 0.5888(0.0024)*** 
c** = -0.1178(0.0061)*** 

P0 = 0.0425(0.0001)*** 

OLp = -0.6079(0.0051)*** 

k 0 = 77.0015(0.1864)*** 
o: K = -9.7247(0.0177)*** 

P0 = 1.4404(0.0035)*** 

OLp = -8.2762(0.0164)*** 

GE on NYSE 

GE on ARCA 

GE on NASDAQ T 

N = 5000, a + = 1.0289 

N = 5000, o+ = 1.0031 

N = 1000*, a + = 1.000000 

P s = 0.534 

Timeout reached. 

Ps = ~ in f 

k 0 = 1.4353(0.0019)*** 

P0 = 0.0432(0.0000)*** 


K 0 = 1.4767(0.0050)*** 

PO = 0.0475(0.0000)*** 


















GE on CBOE 

GE on BATS 

GE on ISE 

N = 5000, a + = 5.2058 

N = 1000*, a + = 1.000000 

N = 1000*, a + = 3.272727 

P s = 0.672 

P s = 1.00 

P s = 0.930 Pm = 0.943 

k 0 = 0.2249(0.0007)*** 

P0 = 0.1067(0.0005)*** 

K 0 = 1.3314(0.0028)*** 

P0 = 0.0343(0.0000)*** 

kq = 32.9135(0.0937)*** 
c* K = -0.2634(0.0164)*** 

P0 = 0.2761(0.0007)*** 
oc p = -0.6564(0.0156)*** 

MKTX on NASDAQ OB 

MKTX on NSE 

MKTX on Chicago 

N = 1303, a+ = 16.015342 

N = 16. a+ = 7.4444 

N = 0, a+ = 0.0000 

Optimization failed 

No data 

No data 

MKTX on NYSE 

MKTX on ARCA 

MKTX on NASDAQ T 

N = 0, a+ = 0.0000 

N = 5000, a+ = 5.0324 

N = 5000, a+ = 3.3373 

No data 

P s = 0.450 P m = 0.403 

P s = 0.457 Pm = 0.442 


k 0 = 0.2160(0.0018)*** 

= 0.0393(0.0007)*** 

= -0.5419(0.0075)*** 

P0 = 0.1022(0.0006)*** 

Pl = 0.1463(0.0019)*** 
oc p = -1.1451(0.0065)*** 

k 0 = 0.0307(0.0008)*** 

K! = 0.0340(0.0012)*** 
oc K = -0.8618(0.0284)*** 

P0 = 0.0475(0.0006)*** 

Pl = 0.0587(0.0018)*** 
otp = -0.7235(0.0261)*** 

MKTX on CBOE 

MKTX on BATS 

MKTX on ISE 

N = 1000, a + = 208.296364 

N = 5000, a + = 14.3535 

N = 5000, a+ = 12.6251 

P s = 0.928 

P s = 0.352 Pm = 0.134 

P s = 0.308 

k 0 = 0.0000(0.0000)*** 

P0 = 0.0000(0.0000)*** 

kq = 1.4282(0.0185)*** 

= 0.0056(0.0001)*** 
o: K = -0.0300(0.0072)*** 

P0 = 1.7276(0.0180)*** 

Pl = 0.0311(0.0006)*** 
ocp = -0.9594(0.0081)*** 

k 0 = 0.0050(0.0000)*** 

P0 = 0.0048(0.0000)*** 

JCOM on NASDAQ OB 

JCOM on NSE 

JCOM on Chicago 

N = 5000, a + = 2.7002 

N = 40, a+ = 9.1136 

N = 0, a+ = 505.0000 

P s = 0.055 Pm = 0.056 

Optimization failed 

No data 

k 0 = 0.5804(0.0379)*** 
ot K = -0.4535(0.0879)*** 

P0 = 1.5632(0.0984)*** 
cc p = -0.4486(0.0870)*** 



JCOM on NYSE 

JCOM on ARCA 

JCOM on NASDAQ T 

IV = 0, a+ = 0.0000 

N = 5000, a + = 1.7622 

N = 5000, a+ = 2.0356 

No data 

P s = 0.169 Pm = 0.177 

P s = 0.376 Pm = 0.382 


k 0 = 9.0255(0.3157)*** 

= 15.3048(1.4759)*** 
ex* = -2.6170(0.0562)*** 

P0 = 7.4636(0.2667)*** 

Pl = 12.9801(1.2344)*** 

OLp = -2.4110(0.0565)*** 

k 0 = 0.4405(0.0081)*** 

«l = 0.3113(0.0067)*** 
q:^ = -1.0299(0.0226)*** 

P0 = 0.3510(0.0053)*** 

Pl = 0.2242(0.0042)*** 
otp = -1.3310(0.0243)*** 
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JCOM on CBOE 

JCOM on BATS 

JCOM on ISE 

N = 5000, a + = 113.4504 

N = 5000, a + = 2.6493 

N = 5000, a+ = 23.2402 

P s = -0.196 Pm = -0.221 

P s = 0.526 Pm = 0.605 

P s = 0.487 P m = 0.390 

k 0 = 0.0007(0.0001)*** 
ki = 0.0876(0.0012)*** 
a K = -1.6064(0.0057)*** 

P0 = 0.0993(0.0170)*** 

Pl = 0.1877(0.0025)*** 
oc p = -0.6117(0.0052)*** 

k 0 = 50.9194(1.6604)*** 

«1 = 25.2161(1.6422)*** 

= -3.1496(0.0218)*** 

P0 = 47.6554(1.5239)*** 

Pl = 22.8284(1.4636)*** 
oc p = -2.6040(0.0205)*** 

k 0 = 0.0167(0.0003)*** 
ki = 0.0095(0.0002)*** 
c* K = -0.4696(0.0063)*** 

P0 = 0.0068(0.0001)*** 

Pl = 0.1036(0.0010)*** 
cc p = -0.7617(0.0034)*** 

ARL on NASDAQ OB 

ARL on NSE 

ARL on Chicago 

N = 0. 3+ = 0.0000 

N = 0, a+ = 0.0000 

N = 0, a+ = 0.0000 

No data 

No data 

No data 

ARL on NYSE 

ARL on ARCA 

ARL on NASDAQ T 

N = 5000, a+ = 7.3767 

N = 1000, a+ = 38.241818 

N = 3639, a+ = 5051.364227 

Optimization failed 

Optimization failed 

P s = 0.083 



k 0 = 0.0000(0.0000)*** 

P0 = 0.0000(0.0000)*** 

ARL on CBOE 

ARL on BATS 

ARL on ISE 

N = 283, a+ = 501008.022436 

N = 97, a+ = 9354.074766 

N = 283, a+ = 72545.182254 

Optimization failed 

Optimization failed 

Optimization failed 


C.2 GZI Model 


XOM on NASDAQ OB 

XOM on NSE 

XOM on Chicago 

N = 5000, a + = 1.928909, q = 1.005368 

N = 1000*, a + - 36.201818, 

q= 1.433895 

N = 46, a+ = 7.4118, q = 2.3594 

P s = 0.225 Pm = 0.194 

P m = 0.704 

Pm = 0.346 

k 0 = 17.2408(1.3449)*** 

= 0.7506(0.0773)*** 
ot K = -2.6031(0.1214)*** 

P0 = 17.9548(1.3943)*** 

Pl = 1.0378(0.1075)*** 
oc p = -2.6008(0.1272)*** 

?7 = 0.9230(0.0107)*** 

k 0 = 0.3677(0.0160)*** 
ot K = -2.8310(0.0348)*** 

P0 = 3.0327(0.0771)*** 
ctp = -2.8919(0.0279)*** 

77 = 0.9845(0.0012)*** 

k 0 = 0.0256(0.0051)*** 

«l = 0.0210(0.0095)* 
at K = -2.2872(0.2175)*** 

P0 = 0.0081(0.0016)*** 

Pl = 0.2698(0.0721)*** 
oc p = -3.2243(0.1371)*** 

77 = 0.9895(0.0118)*** 

XOM on NYSE 

XOM on ARCA 

XOM on NASDAQ T 

N = 5000, a+ = 1.8062, q = 1.9048 

N = 5000, a+ = 1.0658, q = 1.2073 

N = 5000, a+ = 1.0445, q = 1.2832 

P s = 0.168 

P s = 0.838 P m = 0.837 

P s = 0.738 

« 0 = 2.1825(0.0094)*** 

P0 = 0.4998(0.0021)*** 

77 = 0.8044(0.0012)*** 

k 0 = 40.9699(0.6556)*** 

«1 = 44.1512(2.9864)*** 
ct K = -1.7462(0.1082)*** 

P0 = 5.4837(0.0888)*** 

Pl = 7.7849(0.5395)*** 
cc p = -0.7655(0.1062)*** 

77 = 0.7761(0.0029)*** 

kq = 128.2530(2.4351)*** 

P0 = 18.3973(0.3463)*** 

77 = 0.7620(0.0057)*** 
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XOM on CBOE 

XOM on BATS 

XOM on ISE 

N = 662, a+ = 65.5720, q = 1.4527 

N = 5000, a+ = 1.0104, q = 1.1257 

N = 1000*, a + = 17.153636, 

q = 1.333631 

P s = -0.220 Pm = -0.054 

P s = 0.335 P m = 0.313 

P s = 0.453 Pm = 0.439 

« 0 = 0.2861(0.0419)*** 

«l = 0.0676(0.0100)*** 
oc K = -1.6968(0.0410)*** 
p 0 = 4.4734(0.5142)*** 

Pl = 6.9722(0.9634)*** 
cc p = -2.2068(0.0468)*** 

77 = 0.9361(0.0033)*** 

k 0 = 133.1251(2.4196)*** 

= -4.8374(0.1596)*** 

P0 = 16.8156(0.3043)*** 
cc p = -4.5862(0.1486)*** 

77 = 0.7135(0.0038)*** 

k 0 = 2.3030(0.0567)*** 
at K = -2.7609(0.0150)*** 

P0 = 6.1707(0.0873)*** 
cc p = -2.4294(0.0119)*** 

77 = 0.7971(0.0039)*** 

MSFT on NASDAQ OB 

MSFT on NSE 

MSFT on Chicago 

N = 5000, a+ = 1.0331, q = 4.1941 

N = 5000, a+ = 1.0675, q = 1.7971 

N - 1092, a+ = 7.4010, q = 6.0947 

P s = 0.671 Pm = 0.679 

P s = 0.433 Pm = 0.381 

P s = 0.543 P m = 0.490 

«0 = 0.6863(0.0039)*** 

= -0.4756(0.0199)*** 

P0 = 0.0487(0.0003)*** 
u p = -0.6009(0.0192)*** 

77 = 0.8249(0.0006)*** 

« 0 = 0.0099(0.0003)*** 

«l = 0.2286(0.0105)*** 
a: K = -2.6446(0.0807)*** 

P0 = 0.0079(0.0001)*** 

Pl = 0.1003(0.0040)*** 
oc p = -2.3066(0.0740)*** 

77 = 0.9666(0.0011)*** 

k 0 = 0.0020(0.0001)*** 

= 0.0022(0.0000)*** 

= -0.7520(0.0076)*** 

P0 = 0.0028(0.0000)*** 

Pl = 0.0014(0.0000)*** 
up = -1.9361(0.0256)*** 

77 = 0.9214(0.0006)*** 

MSFT on NYSE 

MSFT on ARCA 

MSFT on NASDAQ T 

N = 0, a+ = 0.0000, q = 0.0000 

N = 5000, o+ = 1.0005, q = 2.9264 

N = 0, a+ = 0.0000, q = 0.0000 

No data. 

Timeout reached 

No data 

MSFT on CBOE 

MSFT on BATS 

MSFT on ISE 

N = 296, a+ = 315.9969, q = 1.5671 

N = 1000*, a + = 1.000000, 

N = 1000*, a+ = 1.503636, 


q= 1.637273 

q = 4.573843 

P s = 0.683 Pm = 0.718 

P s = i.oo 

P s = 0.959 Pm = 0.959 

k 0 = 0.0090(0.0012)*** 

= 0.0091(0.0007)*** 
c* K = -0.3375(0.0483)*** 

P0 = 0.0234(0.0022)*** 

Pl = 0.0253(0.0023)*** 
ocp = -9.2445(0.0873)*** 

-q = 0.8491(0.0034)*** 

k 0 = 3.0505(0.0067)*** 

P0 = 0.0407(0.0001)*** 

77 = 0.9397(0.0004)*** 

k 0 = 899.9914(6.3647)*** 
u K = -6.7585(0.0391)*** 

P0 = 17.2797(0.1224)*** 
up = -3.5550(0.0268)*** 

77 = 1.0000(0.0031)*** 

GE on NASDAQ OB 

GE on NSE 

GE on Chicago 

N = 5000, a + = 1.0480, q = 3.5815 

N = 5000, a + = 1.0847, q = 9.0287 

N = 1311, a+ = 2.4913, q = 6.9945 

P s = 0.736 Pm = 0.737 

P s = 0.329 

P s = 0.790 Pm = 0.790 

«0 = 10.3450(0.0536)*** 
a K = -4.1472(0.0266)*** 

P0 = 0.7183(0.0035)*** 
ctp = -4.2231(0.0290)*** 

77 = 0.8351(0.0005)*** 

k 0 = 1.3172(0.0026)*** 

P0 = 0.0543(0.0001)*** 

77 = 0.9654(0.0003)*** 

= 0.4598(0.0024)*** 
c* K = -1.4734(0.0039)*** 

P0 = 0.0184(0.0001)*** 
ctp = -1.7745(0.0059)*** 

77 = 0.7351(0.0009)*** 
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GE on NYSE 

GE on ARCA 

GE on NASDAQ T 

N = 1000*, a+ = 1.009091, 

q = 12.854730 

N = 5000, a + = 1.0024, q = 4.7237 

N = 5000, a+ = 1.0007, q = 6.1216 

P s = 1.000 

Timeout reached 

Timeout reached 

« 0 = 40.2572(0.1865)*** 

PO = 0.6672(0.0030)*** 

77 = 0.9157(0.0006)*** 



GE on CBOE 

GE on BATS 

GE on ISE 

N = 413, a+ = 2.9824, q = 11.6940 

N = 1000*, a+ = 1.000000, 

q = 2.263397 

N = 1000*, a + = 1.456364, 

q = 6.700976 

P s = 0.218 Pm = 0.227 

« 0 = 1.6704(0.0167)*** 

«1 = 4.1414(0.2666)*** 

= -2.1172(0.0500)*** 

P0 = 0.1456(0.0018)*** 
pi = 9.4944(0.5865)*** 
oc p = -3.3851(0.0461)*** 
rj = 0.8420(0.0028)*** 

P s = i.oo 

k 0 = 2.1701(0.0040)*** 

P0 = 0.0373(0.0000)*** 

77 = 0.9770(0.0003)*** 

P s = 0.360 

k 0 = 0.4013(0.0011)*** 

P0 = 0.0018(0.0000)*** 

77 = 0.6216(0.0003)*** 

MKTX on NASDAQ 
OMX BX 

MKTX on NSE 

N = 4, a+ = 6.0000, q = 1.0000 

MKTX on Chicago 

N = 0, a+ = 0.0000, q = 0.0000 

N = 1, a+ = 61.500000, q = 1.000000 

No data 

No data 

Pm = -0.007 



«0 = 0.1442(0.0000)*** 
a K = 0.0000( — nan) 

P0 = 168187.2977(—nan) 
cx.p = —9.6653 (—nan) 
rj = 1.0000 (—nan) 



MKTX on NYSE 

MKTX on ARCA 

MKTX on NASDAQ T 

N = 0, a+ = 0.0000, q = 0.0000 

N = 4379, a + = 5.5327, q = 1.2060 

N = 0, a+ = 0.0000, q = 0.0000 

No data 

P s = 0.314 Pm = 0.311 

No data 


« 0 = 0.0042(0.0001)*** 

«l = 0.0027(0.0001)*** 
o: K = -0.0393(0.0096)*** 

P0 = 0.0042(0.0001)*** 

Pl = 0.0110(0.0001)*** 
ocp = -0.4967(0.0067)*** 

77 = 0.6193(0.0016)*** 


MKTX on CBOE 

MKTX on BATS 

MKTX on ISE 

N = 100, a + = 1049.7568, q = 1.3153 

N = 1623, a + = 18.0118, q = 1.1758 

N = 590, a + = 26.6271, q = 1.2500 

P s = 0.998 

P s = 0.411 

P s = -1.121 

k 0 = 0.0000(0.0000)*** 

P0 = 0.0012(0.0004)*** 

77 = 1.0000(0.0138)*** 

k 0 = 0.0004(0.0000)*** 

P0 = 0.0005(0.0000)*** 

77 = 0.8852(0.0017)*** 

k 0 = 0 . 0010 ( 0 . 0000 )*** 

P0 = 0.0011(0.0000)*** 

77 = 0.3867(0.0031)*** 
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JCOM on NASDAQ 
OMX BX 

N = 1900, a+ = 51.474163, 

q = 1.203507 

P s = 0.904 Pm = 0.772 

« 0 = 0.0050(0.0004)*** 

«l = 0.0014(0.0002)*** 
o: K = -1.8681(0.1083)*** 

P0 = 0.0193(0.0007)*** 

Pl = 0.0256(0.0028)*** 
oc p = -2.2825(0.1025)*** 

77 = 0.9998(0.0001)*** 

JCOM on NSE 

N = 13, a+ = 8.2667, q = 1.0769 

P s = -0.439 Pm = -0.400 

«q = 2.0000(7ia7i) 

= 2.0000(7ia7i) 
oc K = 0 . 0000 ( 7 icm) 

P0 = 2.0000(7ia7i) 

Pl = 1.0000(7icm) 
otp = 0.0000(7ia7i) 

77 = 1.0000(7iGm) 

JCOM on Chicago 

N = 0, a+ = 0.0000, q = 0.0000 

No data 

JCOM on NYSE 

JCOM on ARCA 

JCOM on NASDAQ T 

N = 0, a+ = 0.0000, q = 0.0000 

N = 5000, a+ = 3.0464, q = 1.1506 

N = 0, a+ = 0.0000, q = 0.0000 

No data 

P s = 0.082 

kq = 0.1442(0.0027)*** 

P0 = 0.1998(0.0037)*** 

77 = 0.5970(0.0062)*** 

No data 

JCOM on CBOE 

JCOM on BATS 

JCOM on ISE 

N = 118, a+ = 42.8000, q = 1.1923 

N = 5000, a+ = 15.0924, q = 1.1036 

N = 2224, a+ = 46.0266, q = 1.1379 

P s = 0.295 Pm = 0.291 

P s = 0.690 P m = 0.653 

P s = -0.839 

« 0 = 0.0271(0.0029)*** 

k 0 = 0.0044(0.0002)*** 

k 0 = 0.0000(0.0000)*** 

= -0.4791(0.0330)*** 

= 0.0102(0.0003)*** 

P0 = 0.0000(0.0000)*** 

P0 = 0.0380(0.0045)*** 

o: K = -0.9379(0.0086)*** 

77 = 0.4636(0.0011)*** 

oc p = -0.3331(0.0354)*** 

P0 = 0.0041(0.0001)*** 


77 = 0.9490(0.0051)*** 

Pl = 0.0635(0.0013)*** 
oc p = -1.4285(0.0100)*** 

77 = 0.7080(0.0022)*** 


ARL on NASDAQ OMX 
BX 

ARL on NSE 

ARL on Chicago 

N = 0, a+ = 0.000000, q = 0.000000 

N = 0, a+ = 0.000000, q = 0.000000 

N = 0, a+ = 0.000000, q = 0.000000 

No data 

No data 

No data 

ARL on NYSE 

ARL on ARCA 

ARL on NASDAQ T 

N = 299, a+ = 7.793313, q = 1.310448 

N = 78, a+ = 16204.7558, q = 1.5889 

N = 90, a+ = 17.1500, q = 1.1569 

P s = -7.992 

Optimization failed 

P s = -0.246 

k 0 = 0.0000(0.0000)*** 


k 0 = 0.0000(0.0000)*** 

P0 = 0.0001(0.0000)*** 


P0 = 4.8238(0.4067)*** 

77 = 0.8833(0.0089)*** 


77 = 1.0000(0.0494)*** 
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ARL on CBOE 

ARL on BATS 

ARL on ISE 

N = 0, a+ = 150.0000, q = 1.0000 

N = 7, o+ = 78.8750, q = 1.1250 

N = 6 , a + = 28469.4286, q = 1.8571 

No data 

P s = 0.997 

P s = 0.994 P m = 0.994 


«0 = 0.0219(0.0019)*** 

P0 = 0.1378(0.0279)*** 

77 = 1.0000(0.0789)*** 

KQ = 2.0000(7ia7l) 

Kl = 2.0000(t2O7i) 
ot n = 0.0000(7ioti) 

PO = 2.0000(7iO7i) 

Pl = 1. 0000 ( 71071 ) 

OLp = 0.0000(71O7l) 

77 = 1.0000(nan) 

ARL on NASDAQ OMX 
BX 

ARL on NSE 

N = 0, a+ = 0.000000, q = 0.000000 

ARL on Chicago 

N = 0, a+ = 0.000000, q = 0.000000 

IV = 0, a+ = 0.000000, q = 0.000000 

No data 

No data 

No data 



ARL on NYSE 

ARL on ARCA 

ARL on NASDAQ T 

N = 299, o+ = 7.793313, q = 1.310448 

N = 78, a + = 16204.7558, q = 1.5889 

N = 90, o+ = 17.1500, q = 1.1569 

P s = -7.992 

P s = 0.998 P m = 0.998 

P s = -0.246 

k 0 = 0.0000(0.0000)*** 

P0 = 0.0001(0.0000)*** 

77 = 0.8833(0.0089) *** 

kq = 2.0000(7iem) 
ki = 2.0000(7ia7i) 
at K = 0.0000(7icm) 

P0 = 2.0000(7ia7i) 

Pl = 1.0000(7icm) 

OLp = 0.0000(720.71) 

77 = 1. 0000 ( 7 ia 7 i) 

k 0 = 0.0000(0.0000)*** 

PO = 4.8238(0.4067)*** 

77 = 1.0000(0.0494)*** 

ARL on CBOE 

ARL on BATS 

ARL on ISE 

N =0,a+ = 150.0000, q = 1.0000 

N = 7, a + = 78.8750, q = 1.1250 

N = 6, a + = 28469.4286, q = 1.8571 

No data 

P s = 0.997 

P s = 0.994 Pm = 0.994 


k 0 = 0.0219(0.0019)*** 

PO = 0.1378(0.0279)*** 

77 = 1.0000(0.0789)*** 

kq — 2.0000(7ia7i) 
ki = 2.0000(nan) 
ol k = 0.0000(72a7l) 

P0 = 2.0000(7ia7i) 

Pl = 1.0000(7ia7i) 

OLp = 0.0000(7ia7l) 

77 = 1. 0000 ( 7 io 7 i) 
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C.3 ZI Model with Trades 


XOM on NASDAQ OB 

XOM on NSE 

XOM on Chicago 

N = 5000, a+ = 1.928909 

N = 1000*, a+ = 36.201818 

N = 46, a+ = 7.4118 

P S = 0.183 Pm = 0.170 

P s = 0.649 Pm = 0.695 

Pm = 0.340 

k 0 = 42.3316(4.0871)*** 
ki = 1.3307(0.1302)*** 
c** = -2.8164(0.1107)*** 
p 0 = 44.2052(4.2340)*** 

Pl = 1.8553(0.1773)*** 
oc p = -2.6764(0.1054)*** 

«0 = 0.3855(0.0167)*** 

= -2.7543(0.0340)*** 

P0 = 3.0505(0.0657)*** 
a p = -2.6481(0.0156)*** 

« 0 = 0.0265(0.0048)*** 

«l = 0.0188(0.0079)** 
q:k = -2.2116(0.2067)*** 

P0 = 0.0084(0.0015)*** 

Pl = 0.2206(0.0529)*** 
cc p = -3.0384(0.1210)*** 

XOM on NYSE 

XOM on ARCA 

XOM on NASDAQ T 

N = 5000, a+ = 1.8062 

N = 5000, a + = 1.0658 

N = 5000, a + - 1.0445 

P s = 0.133 

P s = 0.835 P m = 0.836 

P s = 0.739 Pm = 0.739 

k 0 = 3.5635(0.0135)*** 

P0 = 0.9331(0.0030)*** 

« 0 = 122.6212(1.3671)*** 

= 48.9888(3.6449)*** 

= -1.5287(0.0914)*** 

P0 = 17.0316(0.1812)*** 

Pl = 9.6314(0.7186)*** 
oc p = -0.5917(0.0862)*** 

k 0 = 328.9379(5.7143)*** 
oc K = -1.0492(0.2231)*** 

P0 = 47.4101(0.8032)*** 
cc p = -1.0763(0.2234)*** 

XOM on OBOE 

XOM on BATS 

XOM on ISE 

N = 662, a+ = 65.5720 

N = 5000, a+ = 1.0104 

N = 1000*, a + = 17.153636 

P s = -0.335 P m = -0.069 

P s = 0.599 Pm = 0.598 

P s = 0.417 

k 0 = 0.2892(0.0431)*** 

= 0.0698(0.0109)*** 
ot K = -1.6006(0.0429)*** 

P0 = 4.6186(0.5197)*** 

Pl = 6.4110(0.8990)*** 
ot p = -1.9299(0.0365)*** 

K 0 = 4547.3608(195.8280)*** 
ot K = -3.3412(0.3168)*** 

P0 = 584.8992(25.1806)*** 
otp = -2.9657(0.3183)*** 

k 0 = 0.0529(0.0006)*** 

P0 = 0.3351(0.0020)*** 

MSFT on NASDAQ OB 

MSFT on NSE 

MSFT on Chicago 

N = 5000, a+ = 1.0331 

N = 5000, a+ = 1.0675 

N = 1092, a+ = 7.4010 

P s = 0.665 P m = 0.664 

P s = 0.303 Pm = 0.288 

P s = 0.463 Pm = 0.512 

k 0 = 3.2470(0.0128)*** 
ai K = -0.0807(0.0286)** 

P0 = 0.2541(0.0008)*** 
cc p = -0.3001(0.0270)*** 

k 0 = 0.0167(0.0005)*** 
kx = 0.1939(0.0095)*** 

= -2.4091(0.0680)*** 

P0 = 0.0135(0.0001)*** 

Pl = 0.0823(0.0037)*** 
oc p = -2.0423(0.0669)*** 

k 0 = 0.0012(0.0000)*** 
c* K = -0.3006(0.0072)*** 

P0 = 0.0024(0.0000)*** 
otp = -0.9152(0.0049)*** 

MSFT on NYSE 

MSFT on ARCA 

MSFT on NASDAQ T 

N = 0, a+ = 0.0000 

N - 1000*, a+ = 1.000000 

N = 0, a+ = 0.0000 

No data 

Ps = 1.0 

No data 


«0 = 1.4199(0.0037)*** 

P0 = 0.0124(0.0000)*** 
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MSFT on CBOE 

MSFT on BATS 

MSFT on ISE 

N = 296, a+ = 315.9969 

N = 1000, a+ = 1.000000 

N = 1000*, a + = 1.503636 

P s = 0.680 

P S = 1.0 

P s = 0.959 P m = 0.959 

k 0 = 0.0146(0.0004)*** 

PO = 0.0199(0.0006)*** 

K 0 = 1664.8294(16.2923)*** 

P0 = 16.1257(0.1567)*** 

kq = 899.9914(6.1338)*** 

= -6.7585(0.0389)*** 

P0 = 17.2797(0.1171)*** 
oc p = -3.5550(0.0267)*** 

GE on NASDAQ OB 

GE on NSE 

GE on Chicago 

N = 5000, a+ = 1.0480 

N = 5000, a+ = 1.0847 

N = 1311, a+ = 2.4913 

P s = 0.742 

Timeout reached 

P s = 0.790 Pm = 0.790 

k 0 = 68.1041(0.2727)*** 

PO = 5.0341(0.0183)*** 


k 0 = 3.4210(0.0154)*** 
a: K = -2.5542(0.0051)*** 

P0 = 0.1435(0.0006)*** 
ol p = -2.9388(0.0073)*** 

GE on NYSE 

GE on ARCA 

GE on NASDAQ T 

N = 1000* , a+ = 1.009091 

N = 5000, a+ = 1.0024 

N = 5000, o+ = 1.0007 

P s = 1.000 

Timeout reached 

Timeout reached 

«0 = 67.4131(0.2748)*** 

PO = 1.1715(0.0046)*** 



GE on CBOE 

GE on BATS 

GE on ISE 

N = 413, a+ = 2.9824 

N = 1000*, a+ = 1.000000 

N = 1000*, a+ = 1.456364 

P s = 0.222 P m = 0.233 

Ps = 1.0 

P s - 0.417 

k 0 = 1.9455(0.0178)*** 

= 4.4616(0.2844)*** 
at K = -1.8417(0.0465)*** 

PO = 0.1828(0.0017)*** 

Pl = 8.6946(0.5419)*** 
otp = -2.8376(0.0460)*** 

K 0 = 1.4835(0.0031)*** 

P0 = 0.0392(0.0000)*** 

k 0 = 1.1043(0.0026)*** 

P0 = 0.0093(0.0000)*** 

MKTX on NASDAQ OB 

MKTX on NSE 

MKTX on Chicago 

N = 1, a+ = 61.500000 

N = 4, a+ = 6.000000 

N = 0, a+ = 0.0000 

No data 

No data 

No data 

MKTX on NYSE 

MKTX on ARCA 

MKTX on NASDAQ T 

N = 0, a+ = 0.0000 

N = 4379, a + = 5.5327 

N = 0, a+ = 0.0000 

No data 

P s = 0.269 P m = 0.250 

No data 


k 0 = 0.0666(0.0010)*** 

= 0.0167(0.0003)*** 
ex* = -0.4550(0.0080)*** 

P0 = 0.1096(0.0011)*** 

Pl = 0.0708(0.0006)*** 
ol p = -0.8269(0.0048)*** 
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MKTX on CBOE 

MKTX on BATS 

MKTX on ISE 

N = 100, a+ = 1049.7568 

N = 1623, a + = 18.0118 

N = 590, a+ = 26.6271 

P s = 0.996 

P s = 0.394 

P s = -1.942 

k 0 = 0.0000(0.0000) *** 

P0 = 0.0018(0.0002)*** 

= 0.0005(0.0000)*** 

P0 = 0.0012(0.0000)*** 

kq = 8.5760(0.2158)*** 

P0 = 27.1661(0.6893)*** 

JCOM on NASDAQ OB 

JCOM on NSE 

JCOM on Chicago 

N = 1900, a+ = 51.4742 

N = 13, a+ = 8.2667 

N = 0, a+ = 0.0000 

P s = 0.903 Pm = 0.789 

No data 

No data 

k 0 = 0.0054(0.0004)*** 
ot K = -2.5172(0.0804)*** 

P0 = 0.0248(0.0007)*** 
oc p = -2.0220(0.0440)*** 



JCOM on NYSE 

JCOM on ARCA 

JCOM on NASDAQ T 

N = 0. a+ = 0.0000 

N = 5000, a + = 3.0464 

N = 0, a+ = 0.0000 

No data 

P s = 0.024 P m = 0.039 

No data 


kq = 0.9528(0.0247)*** 

= 0.7065(0.0241)*** 
cx K = -0.1111(0.0271)*** 

P0 = 1.5036(0.0341)*** 

Pl = 1.1213(0.0359)*** 

otp = -0.2881(0.0260)*** 


JCOM on CBOE 

JCOM on BATS 

JCOM on ISE 

N = 118, a+ = 42.8000 

N = 5000, a + = 15.0924 

N = 2224, a+ = 46.0266 

P s = 0.270 Pm = 0.251 

P s = 0.662 

P s = -1.287 

k 0 = 0.0250(0.0021)*** 
ex* = -0.3670(0.0301)*** 
p 0 = 0.0362(0.0044)*** 
ot p = -0.1593(0.0367)*** 

k 0 = 0.0027(0.0000)*** 

P0 = 0.0138(0.0001)*** 

KQ = 0.0007(0.0000)*** 

P0 = 0.0048(0.0000)*** 

ARL on NASDAQ OB 

ARL on NSE 

ARL on Chicago 

N = 0, a+ = 0.0000 

N = 0, a+ = 0.0000 

N = 0, a+ = 0.0000 

No data 

No data 

No data 

ARL on NYSE 

ARL on ARCA 

ARL on NASDAQ T 

N = 299, a+ = 7.7933 

N = 78, a+ = 16204.7558 

N = 90, a+ = 17.1500 

P s = -6.946 P m = -6.262 

Optimization failed 

P s = 0.426 P m = 0.489 

k 0 = 0.0001(0.0000)*** 

= 0.0000(0.0000)*** 
a K = -0.1784(0.0485)*** 

P0 = 0.0006(0.0001)*** 

Pl = 0.0006(0.0001)*** 

ccp = -1.1138(0.0814)*** 


kq = 0.0037(0.0014)** 
oi K = -1.7762(0.1290)*** 

P0 = 0.0418(0.0151)** 
oc p = -1.5838(0.1053)*** 
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ARL on CBOE 

ARL on BATS 

ARL on ISE 

N = 0, a+ = 0.0000 

N =7,a+ = 78.8750 

N = 6, a+ = 28469.4286 

No data 

No data 

No data 
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